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Abstract 

We study the single-photon signals expected at LEP in models with a 
very light gravitino. The dominant process is neutralino-gravitino production 
(e + e~ — ► \G) with subsequent neutralino decay via x ~^ jG, giving a 7 + E m i ss 
signal. We first calculate the cross section at arbitrary center-of-mass energies 
and provide new analytic expressions for the differential cross section valid for 
general neutralino compositions. We then consider the constraints on the grav- 
itino mass from LEP 1 and LEP161 single-photon searches, and possible such 
searches at the Tevatron. We show that it is possible to evade the stringent 
LEP 1 limits and still obtain an observable rate at LEP 2, in particular in the 
region of parameter space that may explain the CDF eejj + ET, m iss event. As 
diphoton events from neutralino pair-production would not be kinematically 
accessible in this scenario, the observation of whichever photonic signal will 
discriminate among the various light-gravitino scenarios in the literature. We 
also perform a Monte Carlo simulation of the expected energy and angular dis- 
tributions of the emitted photon, and of the missing invariant mass expected in 
the events. Finally we specialize the results to the case of a recently proposed 
one-parameter no-scale supergravity model. 
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1 Introduction 



Supersymmetric models with a light gravitino (67) have been considered for some 
time 0, |, [|, |[, but interest on them has recently surged || [7| because of their 
ability to explain naturally the puzzling e + e~77 + Ex.miss event observed by the CDF 
Collaboration |§. If the gravitino is the lightest supersymmetric particle (LSP),Q 
the next-to-lightest supersymmetric particle (NLSP, typically the lightest neutralino 
(%), as we will assume here) becomes unstable and eventually decays into a photon 
plus a gravitino (x — > jG) f2|. This decay becomes of experimental interest when 
it happens quickly enough for the photon to be observed in the detector. Because 
the interaction of the gravitino with matter is inversely proportional to the gravitino 
mass, the neutralino lifetime will be short enough for a sufficiently light gravitino: 
tUq < 250 eV ||. On the other hand, the gravitino may not be too light, as otherwise 
it would be copiously produced leading to distinctive signals at colliders that have 
not been observed |3], || or cosmological [10] and astrophysical |TTJ embarrassments: 



ttlq > 10 6 eV. LEP 1 searches strengthen this limit to itlq > 10 3 eV in large regions 
of parameter space, when m x < Mz |~2] . 



Theoretically, light gravitinos are expected in gauge-mediated models of low- 
energy supersymmetry ||, where the gravitino mass is related to the scale of su- 
persymmetry breaking via « 6 x 10~ 5 eV (Asusy/500 GeV) 2 . Special cases of 
gravity-mediated models may also yield light gravitinos, when the scale of local and 
global breaking of supersymmetry are decoupled, as in the context of no-scale su- 
pergravity [7|], in which case ~ (m 1 / 2 /Mp/) p M Ph with the gaugino mass 
scale and p ~ 2 a model-dependent constant. 

Experimental searches for supersymmetry are considerably more sensitive in 
this type of neutralino-unstable supersymmetric models. First of all, the lightest 
observable supersymmetric channel is no longer a pair of charginos but in- 

stead a pair of (the usually lighter) neutralinos (xx)> or if t ne gravitino is light 
enough (m^ < 10~ 3 eV) the neutralino-gravitino channel (xG). These new chan- 
nels allow a deeper exploration into parameter space. Furthermore, because of the 
photonic signature in all supersymmetric processes, it becomes possible to over- 
come the loss of experimental sensitivity that occurs when the daughter leptons 
become too soft (as in chargino pair production when m x ± — m x < 10 GeV or 
m x ± > rrii, > m x ± — 3 GeV) , and therefore absolute lower bounds on sparticle masses 
become experimentally attainable in this class of models. Indeed, diphoton searches 
at LEP161 (i.e., = 161 GeV) |13| have been recently shown [14]] to exclude a 



significant fraction of the parameter space that is preferred by the supersymmetric 
interpretations of the CDF event .0 Ongoing runs at LEP 2 should be able to probe 
even deeper into the remaining preferred region of parameter space. 

1 We assume that R-parity is conserved, as otherwise the decay p — > GK + may occur at an 
unsuppressed rate ||. 

2 We should emphasize that even though we are encouraged by the natural interpretation of the 
CDF event within certain light gravitino scenarios, we believe that such scenarios are interesting in 
their own right and should be fully explored irrespective of the status of the CDF event. 
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Our purpose here is to consider in detail a complementary signal in light- 
gravitino models, namely the associated production of gravitinos with neutralinos.Q 
The resulting single-photon signal has been recently shown to be observable at LEP 2 
in certain range of gravitino masses, but only when the diphoton signal from neu- 
tralino pair production is itself not kinematically accessible |12[| . Therefore, experi- 
mental observation of whichever photonic signal will provide very useful information 
in sorting out the various light-gravitino scenarios in the literature. The gravitino 
mass plays a central role in gravitino-production processes, whose rate is inversely 
proportional to the gravitino mass squared (l/m|,). In contrast, the precise value 
of the gravitino mass plays a minor role in the production of the traditional super- 
symmetric particles, as it determines only the decay length of the neutralino. The 
neutralino-gravitino process of interest at LEP is 

~ ffi / TH? 

cr{e + e- -> X G -> 7 + E miss ) cc — , with 3 = \ 1 * , (1) 

m~ V s 

G 

which provides an experimental handle on the gravitino mass. This process was 
considered originally by Fayet []3| (in the restricted case of a very light photino-like 
neutralino) who noted that the 8 threshold behavior in Eq. ([]]) results from sub- 
tle cancellations among all contributing amplitudes. Dimensional analysis indicates 
that this cross section exceeds electroweak strength when M|/(Mp L m~) > a wea k or 

mg < a w Hl Mf/Mpi ~ lCT 4 eV. In the context of LEP 1, this process was revisited 
in the restricted case of a neutralino with a non- negligible zino component, where the 
resonant Z-exchange diagram dominates ||. 

In this paper we first calculate the cross section for the e + e~ — > xG process 
at arbitrary center-of-mass energies and give new analytic expressions for the corre- 
sponding differential cross section (Sec. §). Next we reassess the constraints on the 
gravitino mass in view of the full LEP 1 data set and imposing the preliminary lim- 
its obtained recently from runs at LEP161 (-y/s = 161 GeV), for general neutralino 
compositions (Sec. |3[). We also comment on the potential of analogous searches at 
the Tevatron. We then perform a Monte Carlo simulation of the production and de- 
cay processes leading to the single-photon signal and obtain energy (E 7 ) and angular 
(cos# 7 ) distributions for representative points in parameter space (Sec. |]), and also 
discuss the missing invariant mass distribution expected in the events. We show that 
one may evade the LEP 1 limits and still obtain observable single-photon signals at 
LEP 2, although only when the diphoton signal from neutralino pair production is 
kinematically inaccessible. Finally we specialize our results to the case of our pro- 
posed one-parameter no-scale supergravity model |7|, [TBI] (Sec. |). Our conclusions 
are summarized in Sec. ^. 

3 Gravitino pair production (e + e~ — > GG) exceeds neutralino-gravitino production (e + e~ — > %G) 
only for gravitino masses that have already been excluded experimentally (mQ < 10~ 6 cV). The 
single-photon signal is further suppressed by a factor of a from the radiated photon. 
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2 The e + e — > xG process 

The Feynman diagrams for neutralino-gravitino associated production at LEP are 
shown in Fig. Ill and include s-channel 7 and Z exchange, and t- and w-channel 
selectron {e^i) exchange. At the Z peak one expects the s-channel Z exchange 
diagram to dominate, and one may simply calculate the amplitude for Z — > \G 
decay ||. This result is accurate as long as the neutralino has a non-negligible 
zino component. However, for photino-like neutralinos the other diagrams become 
important. This is also the case for any neutralino composition for center-of-mass 
energies away from the Z peak (y's > Mz). To deal with all cases at once, we 
perform the complete calculation of all diagrams contributing to e + e~ — > \G- We 
first present the general form of the differential cross section and later specialize the 
result for the particular case of a photino-like neutralino in order to expound on 
certain theoretical issues and generalizations of our results. 

In calculating interactions of gravitinos with matter one can proceed in one of 
two ways. One may calculate with the full couplings in the supergravity Lagrangian, 
in which case the vertices of interest are given by |], [ToJ 

e Gf, oc P R>L f- ; 7ct 7 G M oc — p? [7,, 7/ , , (2) 

where the (spin-| goldstino component of the) gravitino field is G M cx d^ip/rriQ, and 
M = 2.4 x 10 18 GeV is the appropriately scaled Planck mass. Alternatively one may 
calculate using a set of much-simplified effective goldstino couplings || 

2 2 

The full and effective couplings give the same results for the cross sections of pro- 
cesses where the typical bad high-energy behavior of the gravitational amplitudes is 
cancelled completely by the diagrams involving only gravitinos and regular super- 
symmetric particles. This is the case for the e + e~ — > \G process in hand, where we 
have verified (in the pure photino limit) that both ways of doing the calculation give 
identical results. For a derivation and explanation of the meaning of the effective 



couplings see Ref. ||10|| . The simplification of using the effective couplings is not bene- 
ficial in other processes, such as gravitino pair-production, where diagrams including 
graviton exchanges must also be included to cancel the bad high-energy behavior of 



the amplitudes [|3], [K|. Also, it is not clear whether this simplification may be used in 
the case of broken gauge symmetries, and therefore we have used the full couplings 
in the case of neutralino compositions other than pure photino, where the s-channel 
Z-exchange amplitude must be taken into account. 
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2.1 General case 

The differential cross section for a general neutralino composition is given by 

da (s — m 2 ) F(s,t,u) 
dcosO = 32tts 2 6(Mm 6 ) 2 ' 

where as usual we define 



(4) 



t = -I( s - m 2)(i_cos0) , (5) 
u = -|(s-mj)(l + cos0) . (6) 

The function F(s, t, u) receives contributions from each amplitude squared and various 
interference terms (some of which vanish). In an obvious notation, we find 

F = F 77 + F tt + F uu + F Z z + Fjt + F JU + F zt + F Zu + F jZ , (7) 

where 

, , x9 2s(s-m 2 )(t 2 + w 2 ) 
F 77 = (N[ ie ) 2 ^ ^ (8) 

2 ^2 + \/ -A 4.2 1^2 



t 2 (m 2 -t)(-t) t 2 (m 2 -t)(- t) 

(*-<) 2 1 L) (t-mlY 



Fu = {X R f v * ' 9 y +(X L ) 2 * " ' (9) 



2 u 2 (m 2 — u) 2 u 2 (m 2 — u)(—u) 
F uu = [Xr) 2 —z h (X L ) (10) 



F ( N , 9 \ 2 (4 + c 2 ) 2s(s-m 2 x )(t 2 + u 2 ) 

ZZ \ l2 cose w ) 2 (s — M|) 2 + (T Z M Z ) 2 { J 

F lt - ( N neX R ) s{t _ ml) -(N u eX L ) s{t _ ml) (12) 
s{u-ml R ) s(u-mj L ) 



cos^y [( S -M 2 ) 2 + (r z M z ) 2 ](t-mlJ 

fiV' c X g ^ (-t)(2st 2 )(s - M 2 ) 

" 1 ' 12 ' ^cosOwJ [(s - M 2 ) 2 + (r z M z ) 2 ](t - mfj ll4J 



Fz u = -[n' 12 c r X r 



g \ (—■ u)(2sm 2 )(s — M 



cos^y [(s - M 2 ) 2 + (r z M z ) 2 ]( M - m|J 

, v , v 9 \ (-u)(2su 2 )(s - Ml) 

+ ( N 12 c L X L ) — A/f2 , 2 | m u , 21 ,„, — -2-T ( 15 ) 



cos^/ [(s - M 2 ) 2 + (r z M z ) 2 ]( M - 



m 



■UN'N'r 9 ^ ^ + C ^ 2 S ( S -m 2 )(t 2 + M 2 )( S -M 2 ) 
Z "^^^J— 2 s[(s - M 2 ) 2 + (r,M z ) 2 ] (16) 
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In these expressions we have used the following e-e^irX (Xr,l) and e-e-Z {cr,l) 



couplings [16 



X R 
X L 

cr 
cl 



N , ,9_^0w 
u cos6»l 



7 \V 
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-N' n e-N' l2 -^{\-sxn 2 9 w ) , 



cos 9 



if 



sin 9w , 



•~ + sin 2 9 W 



(17) 

(18) 

(19) 
(20) 



where N' u and N' 12 denote the photino and zino components of the neutralino re- 
spectively. Indeed, the lightest neutralino may be written in two equivalent ways 
PI: 



x = N[ 1 j + N[ 2 Z + N 13 H° + N U H^ 
X = N n B + N 12 W 3 + N 13 H° + N 1A H° 2 



related by 

N' n = N n cos 9 W + N 12 sin 9 W 



N[ 2 = -N u sin 9 W + N 12 cos 9 



7 \r 



(21) 
(22) 



(23) 



2.2 Pure photino case 

This special case is useful in order to expose various subtleties in the calculation 
that become less apparent (although they are still present) in the case of a general 
neutralino composition. The e + e~ —>■ 7G case is also important because the result can 
be readily taken over to the case of gluino-gravitino production in quark- ant iquark 
annihilation at hadron colliders (qq — > gG). 

In this special case the couplings of the neutralino (photino: = 1, N[ 2 = 0) 
to matter are very simple: X R = e = —X L , F zz = F zt = F Zu = F lZ = 0, and 



where 



2s(s-m 2 )(t 2 + u 2 ) , 

Fl = — i ^ J - (25) 

' s 

- t 2 (m 2 x -t)(-t) u 2 (m 2 x -u)(-u) t 2 (-2st) u 2 (-2su) 

F R = u *2\2 + ^2 \2 + Tu r^T + TTT. T^T ( 26 ) 



(t — ml R ) 2 (u — m| ) 2 s(t — iTi 2 R ) s(u — m 
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^ = t 2 (m 2 x -t)(-t) | u 2 (m 2 x -u)(-u) | t 2 {-2st) | u\-2su) 
L (t~ m i L ) 2 (u-mjj 2 s(t-mj L ) s(u-mj L ) 

With this relatively simple expression we can verify certain expected behaviors 
of the cross section. First, in the limit of unbroken supersymmetry: m x — > m 7 = 0, 
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rrie L R — > m e ~ 0, s + t + u = m x — > 0, one can readily verify from the above equations 
that — > 0. The vanishing of the cross section in this limit is expected as the spin- 
al component of the gravitino (the goldstino) becomes an unphysical particle when 
supersymmetry is unbroken, as it is no longer absorbed by the gravitino to become 
massive. 

A related manifestation of this phenomenon can be exposed by studying the 
threshold behavior of the cross section. The spin-| (goldstino) component of the 
gravitino is essentially obtained by taking the derivative of the full gravitino field, 
thus making the goldstino couplings be proportional to the goldstino momentum 
(/c M ). At threshold —>■ and there is an additional suppression of the cross section 
besides the kinematical one. Threshold corresponds to the limit s — > m x and therefore 
from Eqs. (|5|,|B|) t,u go to zero as (s — m x )- In the above expression for F^, one can 
see that near threshold each term is proportional to (s — m x ) 3 , which combined with 
the (s — mi) term from the phase space integration [see Eq. (^)] yields a cross section 

proportional to f3 8 with (3 = ^Jl — m x /s. 

The above results were originally obtained by Fayet [|J based on a calculation 
of the cross section using the effective couplings (Eq. @), and have been obtained 
here for the first time using the full couplings (Eq. @). Such equivalent expression 
for the cross section makes more evident some further properties of the results, and 
we thus give it explicitly here too. The expression for F^ using the effective couplings 
becomes 



2s(s — m 2 ) + Auts/m 2 



i 4 

+ rn~ 



+ m~ eL 



K-t)H) , ( 
it - m\ R f 

K -*)(-*) 



mj L ) 2 



u)(-u) 



[m 



m\ R f 
u)(-u) 



(u-m%y 



i 2 2 

+ m x m iR 



-2st) 



+ 



'-2.su) 



s (t- m i R ) s(u-ml R ) 

{-2st) (-2su) 
s(t — m| 



s\u-m% L/ 



f28) 



Despite the seemingly different appearances of Fj s and F^, it can be verified (at 
least numerically) that they give identical results. Using F^ s it is immediately ap- 
parent that the cross section vanishes in the unbroken supersymmetry limit (i.e., 
m x ,mg RL — >• 0), as it should. (Note also that for a massless photino (a case of in- 
terest in the early literature) the s-channel diagram does not contribute.) The j3 8 
threshold behavior is not so apparent this time. One can first note that near thresh- 
old F^ G becomes independent of ra^ R L and depends only on m 2 . A little algebra then 

shows that indeed, near threshold, F^ s oc (s — m 2 ) 3 , and thus the same [3 8 threshold 
behavior results, although this time as a result of a cancellation among all of the 
contributing amplitudes. 

The F^g form is also useful in exhibiting the dependence of the cross section 
on the selectron masses. As is evident from Eq. (^), the cross section increases with 
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increasing selectron masses, eventually saturating for very large values of mg. Thus, 
the decoupling theorem still holds (i.e., large values of the sparticle masses have no 
effect), although its specific implementation here is rather peculiar. 

Before moving on to numerical evaluations of the cross sections, let us note 
that the above expressions for the photino cross section (using either the full or 
effective couplings) can be adapted very easily to describe gluino-gravitino production 
in quark-antiquark collisions at the Tevatron or LHC (qq — > gG). In this case the 
process is mediated by s-channel gluon exchange and t-channel q~L,R exchange. One 
needs to replace the e-e^-x (Xr,l) couplings in Eqs. ( |T7| , |i"8"D by those appropriate 
for q-qR,LrX-> one needs to replace the e-e-7 coupling (eMn Eq. ([24])) by the strong 
coupling (g 2 s ), and one needs to insert the appropriate color factor. Of course the 
integration over parton distribution functions also needs to be implemented. (A 
realistic calculation would also include the gluon-fusion channel, which becomes quite 
relevant at LHC energies.) 



3 Experimental constraints 
3.1 LEP 1 

The single-photon signal (7 + E miss ) has been searched for at LEP 1 by various LEP 
Collaborations |17[| . We estimate an upper bound of 0.1 pb on this cross section. This 



estimate is an amalgamation of individual experimental limits with partial LEP 1 
luminosities (~ 100 pb -1 ) and angular acceptance restrictions ( | cos 6*^, | < 0.7). Note 
that the single-photon background at the Z peak (mostly from e + e~vi"-j) is quite 
significant, as otherwise one would naively expect upper bounds of order 3/C < 
0.03 pb. 

A numerical evaluation of the single-photon cross section at LEP 1 versus 
the neutralino mass for mg = 10 _5 eV is shown in Fig. 0, for different choices of 
neutralino composition ('zino': N[ 2 ~ 1; 'bino': Nu = 1, and 'photino': A r { 1 = 1), 
and where we have assumed the typical result B(x — > 7G) = 1.0 In the photino case 
the Z-exchange amplitude is absent (N[ x = 1 =>- N' l2 = 0) and one must also specify 
the selectron masses which mediate the t- and w-channel diagrams; we have taken 
the representative values mg R = mg L = 75, 150 GeV. Increasing the selectron masses 
further leads to only a small increase in the cross section, e.g., at m x = one finds 
a™ z = 1.48, 2.09, 2.36 pb for m s = 150, 300, 1000 GeV, signalling the reaching of the 



decoupling limit for large selectron masses discussed in Sec. pT2 . 

In Fig. ^ we also show (dotted line 'LNZ') the results for a well-motivated 
one-parameter no-scale supergravity model [7[ 115] , which realizes the light gravitino 
scenario that we study here. In this model the neutralino is mostly gaugino, but has 
a small higgsino component at low values of m x , which disappears with increasing 

4 Note that this implies a non- vanishing (possibly small) photino component of the neutralino, as 
would be required in the 'zino' case discussed above. 
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neutralino masses; the neutralino approaches a pure bino at high neutralino masses. 
The selectron masses also vary (increase) continously with the neutralino mass and 
are not degenerate (i.e., m^ L ~ 1.5m e - fl ~ 2m x ). 

This figure makes apparent the constraint on the gravitino mass that arises 
from LEP 1 searches: in some regions of parameter space one must require tjiq > 10 -3 eV 
if m x < Mz^\ To make this result more evident, in Fig. |3] we display the lower bound 
on the gravitino mass versus the neutralino mass, that results from the imposition of 
our estimated upper bound a^ z < 0.1 pb. (The curves that extend beyond m x = Mz 
result from constraints from LEP 2 data, and are discussed below.) Note the depen- 
dence on the selectron mass in the pure photino case. 



3.2 LEP161 

Recent runs of LEP at higher center-of-mass energies have so far yielded no excess 
of single photons over Standard Model expectations. The latest searches at yfs = 
161 GeV have produced upper limits on the single-photon cross section crj; 61 < 1 pb 
rS|. We have evaluated the single-photon cross sections for the neutralino compo- 



sitions used in Fig. |2] at ^fs = 161 GeV. This time all cases depend on the choice 
of selectron masses. The numerical results are shown in Fig. [|, with the experimen- 
tal upper bound denoted by the dashed line. Note that this line extends only for 
m x > Mz, as for m x < Mz the much stronger limits discussed in the previous sec- 
tion apply. Moreover, for = 10 _5 eV (as used in Fig. |j), LEP 1 limits require 
m x > Mz- As the figure makes evident, for m x < Mz the sensitivity to single-photon 
signals at LEP 2 is not competitive with that at LEP 1. 

As discussed above, the cross sections in Fig. |] increase with increasing selec- 
tron masses (saturating at values somewhat larger than the ones shown), and con- 
versely decrease with decreasing selectron masses. The choice of selectron masses also 
affects the near-threshold behavior of the cross section, with light selectron masses 
"delaying" the onset of the /3 8 threshold dependence (see Fig. [|). Note also that the 
photino, bino, and zino cross sections become comparable above the Z pole, when 
the Z-exchange diagram becomes comparable to the other diagrams. In the case of 
the one-parameter model ('LNZ') a peculiar bump appears. This bump is understood 
in terms of the selectron masses that vary continously with the neutralino mass: at 
low values of m x the selectron masses are light and the cross section approaches the 
light fixed-selectron mass curves ('75'); at larger values of m x the selectron masses are 
large and the cross section approaches (and exceeds) the heavy fixed-selectron mass 
curves ('150'). This example brings to light some of the subtle features that might 
arise in realistic models of low-energy supersymmetry. 

In spite of their apparent weakeness, LEP 161 limits on the single-photon cross 
section are useful in constraining the gravitino mass in a neutralino-mass range in- 
accessible at LEP 1. Indeed, decreasing the gravitino mass in Fig. f| by a factor of 3 

5 For reference, in gauge-mediated models of low-energy supersymmetry, such gravitino masses 
correspond to Asusy ~ 3 TeV. 
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will make the cross sections some ten-times larger. The resulting lower bounds on the 
gravitino mass from LEP161 searches are shown in Fig. |3|. This figure shows that, as 
expected, LEP 1 limits dominate for m x < Mz- However, because of the /3 8 threshold 
behavior at y/s = Mz, LEP161 limits 'take over' for neutralino masses slightly below 
Mz, and in the 'photino' case, considerably below M z . 



3.3 Single photons versus diphotons 

It has been made apparent in Fig. |] that for m x < Mz the gravitino mass is con- 
strained to rriQ ^> 10~ 5 eV. If this would indeed be an absolute requirement on the 
gravitino mass [i.e., for all values of m x ) then the cross section for neutralino-gravitino 
production at LEP 2 would be highly suppressed: Fig. | shows that a* 61 < 1 pb 
for itIq = 10~ 5 eV, and a 7 oc m^f '. In other words, if the minimum observable 
single-photon cross section at LEP 2 is ~ 0.1 pb (i.e., for £ ~ 100 pb -1 ), then 
m G ~ 3 x 10~ 5 eV appears to be the limit of sensitivity of LEP 2. 

On the other hand, the process e + e~ — * XX 77 + E m i SS is sensitive to m x < 
\\fs and is independent of the gravitino mass. In light of the single-photon constraints 
on the gravitino mass obtained above, the diphoton process may be observable at 
LEP 2 (i.e., m x < |-y/s<Mz) only if ttlq ^> 10~ 5 eV, and therefore single photons 
will not be simultaneously observable at LEP 2. Conversely, single photons may be 
observable at LEP 2 only if m x > Mz, in which case diphotons may not be observed 
simultaneously (as they require y/s > 2Mz ~ 190 GeV). This dichotomy between 



single-photon and diphoton signals at LEP was first presented in Ref. |12| 



3.4 Other limits 

The above lower limits on itlq are rather significant, and improve considerably on 
previous limits from collider experiments |], |^, |17]] and astrophysical considera- 



tions pll ) as long as m x < Mz- There has also been a recent reassessment |]T9| of the 
hadron collider limits obtained via associated gluino-gravitino production (pp — > gG), 
and via indirect gluino pair-production (pp — > gg) where in addition to the usual su- 
persymmetry QCD diagrams the gravitino is exchanged in the t and u channels. The 
multijet signature of these processes has been contrasted with experimental limits 
from the most recent Tevatron run to show that if gluino pair production is accessible 
at the Tevatron (i.e., rrig < 200 GeV) then a lower limit of itlq > 3 x 10~ 3 eV re- 
sults. This limit and our limits in Fig. |3| may be compared by relating the gluino and 
neutralino masses, as occurs in supergravity models with universal gaugino masses 
at the unification scale: nig ~ 3m x ± ~ 6m x . Therefore, m^<200GeV (as required 



for the bound in Ref. [|19[ to apply) corresponds to m x < 35 GeV. Consulting Fig. y, 
we see that the Tevatron limits are stronger in this neutralino mass range. How- 
ever, the LEP 1(161) limit extends up to m x <Mz(2Mw), which corresponds to 
nig < 550 (960) GeV, which is far from the direct reach of the Tevatron. 

By considering the further processes pp — > gS, gP, where S and P are very light 
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scalar and pseudo-scalar particles associated with the gravitino, the lower bound on 
the gravitino mass becomes much less dependent on the gluino mass and can be taken 
to be ttIq > 3 x 10 _4 eV |[L9| . This lower bound is comparable with those obtained 
above by considering LEP 1 data. However, this result assumes the existence of 
additional light and strongly interacting particles (S and P), an assumption that 
depends on the detailed nature of the mechanism that leads to a very light gravitino. 

3.5 New channels 

Another set of channels of interest at the Tevatron consist of the associated production 
of gravitinos with neutralinos or charginos 

PP^XG^G, (29) 

which have the advantage over pp — > gG of much less phase space suppression. The 
most basic channel is qq — > xG, which leads to a 7 + ET, m iss signal. The cross 
section for this process can be readily obtained from the expressions given in Sec. ^] 
by replacing the initial state electron-positron pairs by quark- ant iquark pairs, the 
exchanged selectrons by squarks, and by integrating the resulting expression over 
parton distribution functions. We have estimated this cross section and find that 
it may be quite significant: up to 85, 25, 15 pb for m x = 50, 75, 100 GeV and = 
10~ 5 eV, in favorable regions of parameter space. In the best case scenario of a 
Tevatron upper limit of 0.1 pb (i.e., 10 events in C = 100 pb -1 ), one may conclude 
that ttiq > (3, 1.6, 1.2) x 10 _4 eV for m x = 50, 75, 100 GeV. Taken at face value, these 
limits are quite competitive with those obtained in Ref. At the moment there 
are no single-photon limits available from CDF nor DO. 

To improve the visibility of the signal, one may want to consider the qq — > x ± G 
channel which, depending on the chargino decay channel, may lead to £ ± + 7 + E Tjmiss 
or 2j + 7 + E^miss signals. The leptonic signal appears particularly promising. For 
all these processes there are some important instrumental backgrounds that need to 
be overcome. For instance pp — *> W — > ez/ e , where the electron is misidentified as 
a photon (i.e., because of limitations in tracking efficiency), leads to a very large 
"single-photon" signal o ■ B(W — > ez/ e ) « 2.4 x 10 3 pb PU| , which may be reduced 
significantly by optimizing the tracking efficiency and making suitable kinematical 
cuts. The other channels mentioned above face similar, although perhaps less severe, 
instrumental backgrounds (e.g., WW — > e + "7" + ET, m iss)- 



4 The single-photon signal 

The total cross section for neutralino-gravitino production has been displayed in Fig. [| 
for a specific center-of-mass energy (^ = 161 GeV) and for some illustrative choices 



of parameter values. The analytic expressions given in Sec. |2.1| allow one to calculate 



these cross sections for arbitrary values of the parameters. In this section we would 
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like to explore some characteristics of the actual signal, i.e., the energy and angular 
distributions of the observable photon and the missing invariant mass distribution in 
the events. 



4.1 Monte Carlo technique 

Our simulation proceeds in a standard way, making use of a 'home-made' Monte 
Carlo event generator. We start in the rest frame of the decaying neutralino, where 
we generate j+G events that are isotropic in this reference frame. Energy-momentum 
conservation requires E' = \p' 7 \ = 2%, which leaves two components of the photon 
momentum to be generated at random (i.e., pL). We then boost the photon momen- 
tum back to the laboratory frame using the neutralino 4-momentum (E x ,p x ), whose 
components are constrained by the kinematics of the e + e~ — > \G process: 

Ex= ^ + I± > | p - x | = ^£_4. m 
x 2 2y/s' IFxl 2 2v/i v ; 

Here we have two components of the neutralino momentum unconstrained (cos 9 X , <p x ). 
For fixed values of these angles we obtain E~ and cosd^ distributions, which are 
purely kinematical effects. The observable distributions are obtained by varying 
(cos# x ,0 x ) and weighing these kinematical distributions with the corresponding dy- 
namical i dctse f ac tors calculable from the expressions given in Sec. 0. 

In what follows we focus on the case of LEP190 (a/s = 190 GeV). First we 
display in Fig. |5| the total cross sections for single-photon production at this center- 
of-mass energy. These should give an idea of the reach in neutralino masses that may 
be accessible at LEP190. As we expect neutralino-gravitino production to be allowed 
only for m x > Mz (to avoid the stringent LEP 1 lower limits on ttiq), we concentrate 
on the following three neutralino mass choices: m x = 100, 125, 150 GeV. To gain some 
insight into the final distributions, we start by displaying the normalized neutralino 
angular distributions - , dcr a as a function of cos# v for m» = 75, 150, 300 GeV, and 

ad cos X a e 

bino (Fig. ||), zino (Fig. [?]), and photino (Fig. H) neutralino compositions. The total 
cross sections for each of the curves can be read off Fig. [5] and for convenience have 
been tabulated in Table [p. As can be seen from Figs. ^, |7|, and || the angular 
distribution varies quite a bit with neutralino mass, although mostly for light selectron 
masses. Note that the angular distributions always remain finite, and generally show 
a preference for the central region. 



4.2 Energy and angular distributions 

The observable photonic energy and angular distributions can be quite unwieldly once 
we allow for the many choices of parameters that we have considered above. An exam- 
ination of all parameter combinations shows that both energy and angular distribu- 
tions are largely insensitive to the neutralino composition, being much more sensitive 
to the mass parameters (i.e., m x ,mi). This result is perhaps not surprising as the 
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Table 1: Total cross sections corresponding to the differential cross sections shown in 
Figs. 6,7,8 at LEP190. All masses in GeV, all cross sections in pb. 



Composition 


m x 


rrie = 75 


rrig = 150 


m g = 300 


bino 


100 


0.34 


0.54 


1.61 




125 


0.32 


0.19 


0.60 




150 


0.20 


0.04 


0.11 


zino 


100 


0.19 


0.49 


1.12 




125 


0.10 


0.17 


0.42 




150 


0.06 


0.03 


0.08 


photino 


100 


0.37 


0.52 


1.64 




125 


0.37 


0.18 


0.61 




150 


0.23 


0.04 


0.12 



observable distributions of relativistic particles are dominated by kinematical effects 
which depend crucially on the mass parameters. Thus, for brevity we show only the re- 
sult in the bino case which, in any event, is representative of typical supergravity mod- 
els. The energy (E y ) and angular (cos# 7 ) distributions for m x = 100, 125, 150 GeV 
are shown in Fig. |9] and Fig. |H] respectively, for mg = 75, 150, 300 GeV. 

The energy distributions (Fig. |9|) show a significant m x and mg dependence. 
As the neutralino mass grows, it tends to produce harder photons. In fact, it is not 
hard to show that in the decay x ~^ lG, with a neutralino energy and momentum as 
in Eq. (]30|) , the photon energy is restricted to the interval 

ml J~s , x 

27s <K ' K 2- (31 > 

as faithfully reproduced in the simulations. (Near threshold (m x « y^s) the photon 
carries away half of the center-of-mass energy.) These distributions show that any 
given single-photon candidate energy (E 7 ) implies an upper bound on the possible 
neutralino masses consistent with the candidate event, 



m 



x 



< J2^~s £ 7 . (32) 



The photonic angular distributions (Fig. |lOf ) are peaked in the forward and 
backward directions, even more so as the neutralino becomes heavier. The selectron 
mass has an interesting effect. In the case of m x = 100 GeV, from Fig. |6] we see 
that for heavy selectron masses the neutralino angular distribution is fairly flat, and 
therefore the photonic distributions should reflect only kinematical effects, as they do 
(i.e., peaked in the forward and backward directions). For light selectron masses, the 
neutralino avoids the forward and backward directions, and the kinematical effect on 
the photons is diminished. 
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4.3 Missing mass distribution 

The dominant background to the neutralino-gravitino signal is single radiative return 
to the Z: e + e~ — > 7Z — > 7^, where the photon is radiated off the initial state and 
the Z boson tends to be on shell. The most distinctive signature for this background 
process appears in the missing invariant mass M m ; ss = yj -E miss — p m i SS distribution, 
which is strongly peaked at M miss w M z . 

The missing mass distribution for the signal can be easily determined, as in 
this case the missing energy and missing momentum are given by 

E m i ss = y/s 5 p m j ss = I p~f \ = Ey , (33) 

and therefore 

M miss = ^s{^~s-2E 1 ) . (34) 

The allowed range of M miss is obtained by inserting the range of photonic energies in 
Eq. (|3~ID ; we obtain 

< Mmiss < yjs-m* . (35) 

Histograms showing missing mass distributions fall in the range specified by Eq. (|35|) . 
and otherwise favor the upper end of the M m j SS range (corresponding to the lower end 
of the photonic energy range). For brevity, we display these distributions only in the 
one-parameter model example discussed in Sec. |5| below. 

We note in passing that the complementary diphoton events from neutralino 
pair-production have a background (e + e~ — > 77.Z — > ^juu) that also peaks for 
M m i SS ~ Mz @. In this case the missing mass in the diphoton signal varies from 
zero up to a maximum value of \{\fs + y 1 s — 4m^), in contrast with the result for 
the single-photon signal in Eq. (^). The M miss distributions of the single-photon and 
diphoton signals differ not only in range, but also in shape. 

5 One- parameter model example 

It should have become clear from the discussion in Sec. |], that the signals to be 
searched for experimentally can have a wide range of characteristics because of the 
variations in the underlying parameters describing the neutralino-gravitino process. 
In reality, the model of supersymmetry that describes nature will have all its mass 
parameters correlated in some way, and the actual observations may be a bewildering 
composite of the many curves shown above. To exemplify this situation, in this section 
we specialize our results to the case of the one-parameter no-scale supergravity model 
that has been mentioned at various places in the preceding discussion. We have 
already shown the corresponding single-photon cross sections at LEP 1, LEP161, 
and LEP 190, as the dashed lines in Figs. 0,f|,|5] respectively. The cross sections for 
\fs > M z show a peculiar bump that, as discussed in Sec. |S7^, can be traced back to 
the fact that the selectron masses vary with the neutralino mass. 
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As a first step towards obtaining the angular and energy photonic distributions, 
in Fig. O we show the normalized neutralino angular distributions - , da „ for \/s = 

° I I ° a a cos X V 

190 GeV and m x = 60, 80, 100, 120, 140 GeV. The total cross sections in each of these 
cases are a = 1.2, 1.3, 1.0, 0.6, 0.2 pb. Note how relatively flat the angular distributions 
are: no more than a 10% variation. This is to be contrasted with the wide range of 
variability observed in the generic cases shown in Figs. || [7|, and ||. In fact, the 
results in the one-parameter model resemble those in the generic models when the 
selectron mass is large (i.e., the dotted lines in those figures, which correspond to 
rrie = 300 GeV). This is to be expected as in the one-parameter model one has 
mg L ~ 1.5 mg R ~ 2m x , indicating increasingly heavier selectrons. 

Following the method outlined in Sec. |], in Figs. [12] and [13] we display the pho- 
tonic energy and angular distributions at LEP190 for three representative neutralino 
masses (m x = 100, 120, 140 GeV). The energy distributions show the same restrictive 
photon energy behavior as predicted by Eq. ([H]). The angular distributions are also 
peaked in the forward and backward directions. 

Finally we consider the missing mass distributions, which are obtained from 
Eq. fl3Hp, and are shown in Fig. O. We note the range of M miss , as prescribed by 



Eq. (^), and the tendency to favor missing mass values toward the upper end of 
the allowed interval. For the neutralino mass choices shown, the missing mass shows 
a distinct preference to be larger than Mz- (This is in contrast with the M n 
distribution in diphoton events, which is more evenly distributed.) 



1 miss 



6 Conclusions 

In this paper we have attempted to study in some detail the physics of supersymmetric 
models with a gravitino light enough that it can be produced directly at collider ex- 
periments. Our discussion has centered mainly around LEP, from where the strongest 
constraints can be obtained at the moment. We have nonetheless outlined the corre- 
sponding program to be followed at the Tevatron, where instrumental backgrounds 
make identification of the single-photon signal a more challenging task. 

We have provided new and explicitly analytical expressions for neutralino- 
gravitino differential cross sections at e + e~ colliders and have discussed some of the 
theoretical subtleties involved in the calculation and some of the peculiar parameter 
dependences of the cross section. We have used our expressions to obtain new lower 
bounds from LEP 1 data on the gravitino mass for m x < Mz- Weaker limits from 
LEP 2 are obtained at higher neutralino masses. Our study includes a Monte Carlo 
simulation of the single-photon signal, which should be helpful in the experimental 
analyses that are just now getting underway. We have also specialized the results to 
our one-parameter no-scale supergravity model, where the signals can be analyzed 
much more simply because of the tight correlations between the model parameters. 
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Figure 1: Feynman diagrams for neutralino-gravitino production at LEP: (a) s- 
channel 7 and Z exchange, and (b) t-channel selectron exchange (e^x,). Additional 
M-channel diagram not shown. 
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Figure 2: Single-photon cross sections (in pb) from neutralino-gravitino production 
at LEP 1 versus the neutralino mass (m x ) for mg = 10 -5 eV and various neutralino 
compositions. The 'photino' curves depend on the selectron mass (75,150). The 
cross sections scale like a oc m~ 2 - The dashed line represents the estimated LEP 1 
upper limit. Also shown is the result for a one-parameter no-scale supergravity model 
('LNZ'). 
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Figure 3: Lower bounds on the gravitino mass (in eV) as a function of the neutralino 
mass (m x ) that result from single-photon searches (7 + E m i SS ) at LEP 1 and LEP161. 
In the 'photino' case at LEP 1 and the 'photino', 'zino', and 'bino' cases at LEP161, 
the selectron mass influences the results. We have chosen m g = 75, 150, 300 GeV, 
denoted by solid, dashed, and dotted lines respectively. Also shown are the bounds 
in a one-parameter no-scale supergravity model ('LNZ'). 
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Figure 4: Single-photon cross sections (in pb) from neutralino-gravitino production at 
LEP 161 versus the neutralino mass (m x ) for mg = 10 -5 eV and various neutralino 
compositions. The solid curves have a fixed value for the selectron mass (75,150), 
whereas the dotted curve corresponds to a one-parameter no-scale supergravity model 
where the selectron masses vary continously with the neutralino mass. The cross 
sections scale like a oc m~ 2 . The preliminary LEP161 upper limit is indicated. 
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Figure 5: Single-photon cross sections (in pb) from neutralino-gravitino production at 
LEP 190 versus the neutralino mass (m x ) for mg = 10 _5 eV and various neutralino 
compositions. The solid curves have a fixed value for the selectron mass (75,150), 
whereas the dotted curve corresponds to a one-parameter no-scale supergravity model 
where the selectron masses vary continously with the neutralino mass. The cross 
sections scale like a oc m~ 2 . 
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Figure 6: Normalized angular distribution of neutralinos of bino composition in 
neutralino-gravitino production at LEP190 for m x = 100, 125, 150 GeV and wig = 
75 (solid), 150 (dashed), 300 (dots) GeV. 
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m x =100 GeV, Vs=190 GeV (zino) 




Figure 7: Normalized angular distribution of neutralinos of zino composition in 
neutralino-gravitino production at LEP190 for m x = 100, 125, 150 GeV and wig = 
75 (solid), 150 (dashed), 300 (dots) GeV. 
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Figure 8: Normalized angular distribution of neutralinos of photino composition in 
neutralino-gravitino production at LEP190 for m x = 100, 125, 150 GeV and wig = 
75 (solid), 150 (dashed), 300 (dots) GeV. 
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Figure 9: Photonic energy distributions in neutralino(bino)-gravitino production at 
LEP190 for m x = 100, 125, 150 GeV and m- e = 75 (solid), 150 (dashed), and 300 
(dots) GeV. 
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Figure 10: Photonic angular distributions in neutralino(bino)-gravitino production 
at LEP190 for m x = 100, 125, 150 GeV and m~ e = 75 (solid), 150 (dashed), and 300 
(dots) GeV. 
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Figure 11: Normalized angular distribution of neutralinos in neutralino-gravitino pro- 
duction at LEP190 for m x = 60, 80, 100, 120, 140 GeV in a one-parameter no-scale 
super gravity model. 
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Figure 12: Photonic energy distributions at LEP190 in a one-parameter no-scale 
supergravity model for m x = 100, 120, 140 GeV. 
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Figure 13: Photonic angular distributions at LEP190 in a one-parameter no-scale 
supergravity model for m x = 100, 120, 140 GeV. 
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Figure 14: Missing invariant mass distributions at LEP190 in a one-parameter no- 
scale supergravity model for m x = 100, 120, 140 GeV. 
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